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Abstract 

The detection of GRB 060218 at z=0.033 by Swift within 1.5 years of operation, 
O ■ together with the detection of GRB 980425 at z = 0.0085 by BeppoSAX, suggest that 

i these low-luminosity (LL) GRBs have a much higher event rate than the canonical 

high-luminosity (HL) GRBs, and they form a distinct new component in the GRB 
luminosity function. We explore the contribution of this previously neglected GRB 
population to the diffuse neutrino background within the internal shock model and 
I compare it with that of the canonical HL population. By considering a wide range of 

' distributions of various parameters (e.g. luminosity, spectral break energy, duration, 

^ . variability time, Lorentz factor, redshift) for both populations, we find that although 

Tlj" I it is difficult to detect neutrinos from the individual LL GRBs, the contribution of 

■ the LL population to the diffuse neutrino background is more than the HL population 
\^ . above about 10^ GeV. 

o : 

^ ■ PACS numbers: 95.85.Ry, 98.70.Rz, 98.70. Sa 

. Keywords: Neutrinos, Gamma Ray Bursts. 

P : 1 Introduction 

+3 . 

■ The possibiUty of high energy neutrino emission from GRBs in the burst and the afterglow 
^ ■ phases has been studied earlier by different groups [1, 2, 3, 4]. It has been found that only 

^ . the nearby high luminosity (HL) GRBs (e.g. GRB 030329) are expected to be detected 

individually by neutrino telescopes like ICECUBE [3, 5]. In general the diffuse neutrino 
background from the whole GRB population is of great observational interests. Recently a 
low-luminosity (LL), long-duration event, GRB 060218/SN 2006aj was detected by Swift [6] 
at a redshift of 0.0331. Earlier, another low luminosity event, GRB 980425/SN 1998bw was 
detected by BeppoSAX at an even lower redshift of 0.0085 [7, 8]. These two LL events share 
some common characteristics such as low luminosity, long duration, and low isotropic 7-ray 
energy. More importantly, their detections at very low redshifts within a relatively short 
period of time (1.5 years for Swift and 6 years for BeppoSAX) imply that they have a much 
higher event rate than the canonical HL GRBs [9, 10, 11, 12]. A straightforward estimate of 
the local rate of LL-GRBs p^^ can be derived from 

ryBeppo r^Swift 

Po^V{z < 0.033) ( r^^PP° + T^^ift\ _ 2, (1) 



where V{z < 0.033) ~ 1.2 x lO^^Gpc^ is the volume enclosed by z = 0.033, Vl^^pp" = 0.123 
and Q^'^^f^ = 1.33 are the solid angles of the GRBM on board BeppoSAX and the BAT on 
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board Swift respectively. The durations of observation are T^^^° ~ 6 yr and T^'^'^f* ~ 1-5 yr 
(when this paper has been written) for BeppoSAX and Swift, respectively. The number 2 on 
the right hand side of Eq.(l) accounts for the detection of GRB 980425 and GRB 060218. 
This gives a rough estimate of p^^ as 800Gpc~^yr~^ Several groups independently estimated 
this event rate[9, 11], e.g. ~ (300 — 1000)Gpc"'Vr~^, which are comparable to the rate of 
LL GRBs 550tI^gGpc-Vr"^ obtained by Liang et al.[12]. These are all much higher than 
the event rate of HL GRBs, pf^ ~ lGpc~^yr~^[13], and its simple extrapolation to the LL 
regime assuming a same population, e.g. Pq ~ 10Gpc~^yr~^[14]. If one assumes that the 
two LL GRBs within z < 0.033 belong to the population of HL GRBs, then the expected 
number of detection of GRBs within z < 0.033 in the BeppoSAX and Swift era is nearly 
0.0026 for ~ lGpc-V~S and is 0.026 for po ~ lOGpc-^r^^ The Poisson probability 
of detecting two events during the observational period is 3.4 x 10~^ for the former and 
3.3 X 10~^ for the latter. Hence, the detection of GRB 060218 strongly suggests the presence 
of a distinct population (LL), which is not a subset of the canonical HL population. This 
implies two distinct components in the GRB luminosity function[12]. One may characterize 
the luminosity function of each population by a broken power law, i.e. 



The spectral indices and break luminosities of this function for LL and HL populations have 
been derived/constrained in [12]. The values of cti, a2, are 0.6, 4.5, 10''^erg/sec and 1.15, 
2.5, lO^^erg/sec for LL and HL GRBs respectively. The values of 5*0 depend on the upper 
and lower limits in the values of luminosity for each population. 

Although the high energy neutrino fiux of these individual LL GRBs is low, the contribution 
to the diffuse neutrino ffux by the whole LL population is non-negligible due to their very 
high event rate. In this paper we explore the neutrino emission of this hitherto ignored 
LL-GRB population, and compare the contributions of the LL- and HL- GRB populations 
to the diffuse neutrino background. 



2 Diffuse Neutrino Spectra from HL and LL GRBs 



We consider the prompt high energy neutrino emission from GRBs in the burst phase due to 
internal shocks. Protons and photons interact to produce pions which subsequently decay to 
neutrinos, p+7 — > tt"^ — > /x+ + z/^ — > e^ + z/e + ^^' + '^M- This is the dominant mechanism of high 
energy neutrino production inside GRBs. We derive the neutrino spectrum using the power 
law nature of photon spectrum as observed by different satellite experiments BATSE [16], 
HETE [15], Swift [17]. GRB 060218 is an X-ray ffash (XRF)[6]. Sakamoto et al.[18] have 
discussed in detail about the characteristics of 45 GRBs observed by HETE-2 (High Energy 
Transient Explorer) [15], some of them were XRFs. They found no statistically significant 
deviation in the distribution of photon spectral indices depending on the type of the source. 
We therefore describe the comoving-frame photon spectrum for both HL and LL GRBs as 




(2) 




(3) 
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The break energy of the photon spectrum in the comoving frame ebr,c is related to its value in 
the somxc rest frame as e^^ = ^^br,c, where F is the bulk Lorentz factor. The normalisation 
constant A is related to the internal energy density U by, 

7- r 71-2 

^ = 7^3^ (4) 

1-72-2 7i-2i 

where 71 < 2, and 72 > 2. The threshold energy for pair production (7 + 7 — > e+ + e~) 
in the comoving frame is ~ 1 MeV. Since the pair production cross section is much larger 
than p7 interaction cross section, the upper cut-off energy in the photon spectrum which 
produces high energy neutrinos should be limited by the threshold energy of pair production. 
Photon energy 1 MeV in the comoving frame corresponds to 300 MeV and 10 MeV in the 
source rest frame for HL and LL GRBs, with Lorentz factor of the order of 300 and 10, 
respectively. For HL GRBs, 300 MeV photons can produce neutrinos with energy 3000 GeV 
in the source rest frame in py interactions through A resonance. For LL GRBs, 10 MeV 
photons can produce neutrinos with energy 100 GeV in the source rest frame. The break 
energies of LL and HL GRBs are assumed to be of the order of keV and MeV, respectively, 
in the source rest frame. Thus the upper limits in photon energies are much higher than 
their break energies for both LL and HL GRBs. The average comoving time is longer 
than the average synchrotron cooling time of protons for HL GRBs, and for LL GRBs they 
are comparable. As we discuss later in detail the average maximum neutrino energies are 
generally constrained by synchrotron cooling of protons. The corresponding lower limits 
in photon energies for p7 interactions are much lower than photon spectral break energies 
for both HL and LL populations. The expression for parameter A has been derived in 
Eq.(4) using the fact that the minimum photon energies are much lower than the break 
energies and the maximum photon energies are much higher than the break energies. The 
threshold energies for synchrotron self absorption [19] of photons are of the order of eV for 
HL GRBs and are much lower for LL GRBs. The lower limit in photon energies derived by 
the maximum proton energies (applying the synchrotron cooling constraint) is higher than 
the typical synchrotron self absorption energies for LL GRBs, and this is also generally the 
case for HL GRBs. If we take 71 = 1. 72 = 2 and a photon energy range of 30KeV to 3MeV 
in the source rest frame, then we get A = as derived in [1]. 

In the comoving frame, the energy loss time (t^r) of a proton of energy Ep^^ due to pion 
production can be derived from 

The variability time is of the order of milliseconds for HL GRBs. For LL GRBs, observations 
show that their lightcurves are generally very smooth; for GRB 060218 the width of the broad 
peak was 431 sec [6]. In p7 interactions the maximum contribution to neutrino production 
comes from the A resonance as discussed in [1]. For 71 < 2 and 72 > 2, we obtain the 
fractional energy loss of protons to pions as a function of proton energy Ep [= E°^'^{1 + z)] 
in the source rest frame. 



1.34T2- 



E„ < E„ 



U{Ep) — fo \ 1.34I1-1 (6) 
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where, f^^ = ^4 r^''"' r 1 ^ i 1 for HL GRBs and f^^ = o-J29L^,47 1 for LL 

GRBs. Here ^^,47 = L^/(10^^ ergs s~^) and i^,2 = i^/(10^ s). For, 71 = 1 and 72 = 2 we get 
an expression of fTr{Ep) for photon energy range 30keV to 3MeV similar to that derived in 

[1], 

f _ n o -^7.51 J Ep/ Epb Ep < Epb . s 

300''f,-3f6r,MeV [ -> J^pb 

where Fsoo = r/300, Fio = F/IO denote typical Lorentz factors of HL and LL GRBs, 
respectively [20, 21]; eb^.MeV and ehr,kev denote the break energy e^^ in units of MeV and 
keV, respectively. The choice of a low e^r for LL GRBs is based on the fact that GRB 
060218 has e^r ~ 5 keV [6], which well satisfy the so-called Amati- relation [22, 23] 



Shr , . / Elf° X 0.51±0.01 , , 

^ = (3.64 ± 0.04 ( ^— ) (8 

lOOfcey ^ ^^7.9 X lO^Ser^^* ^ ^ 

Another LL GRB 980425 apparently does not satisfy this correlation. However, Ghisellini 
et al. [24] argue that GRB 980425 might not be an outlier. In particular, they suspect 
that BeppoSAX only recorded a short duration of hard emission, and it might be similar to 
GRB 060218 should it be detected by an instrument similar to Swift that could record soft 
emission simultaneously. We therefore assume the Amati-relation for all bursts in this paper 
[25], but caution that the diffuse neutrino flux of the LL-component would be overestimated 
(say, a factor of 2) if a fraction (say, 50%) of LL-GRBs do not satisfy the Amati relation. It is 
worth commenting that if most LL-GRBs satisfy the Amati-relation, they are not detectable 
by BATSE, as is the case of GRB 060218 [26], so that they do not contribute to the BATSE 
population. The proton break energy corresponding to the break energy in the photon 
spectrum can be expressed as 

Ef^ = 1.3 X 10'rl,^{ebr,Mev)~'GeY. (9) 

for HL GRBs or E^^ = 1.45 x WTloi€br,kev)~^ GeV for LL GRBs. In p7 interactions both 
7r° and 7r+ can be produced with equal probabilities. 7r+ gets on the average 20% of the 
proton energy and if the final state leptons share the pion energy equally then each neutrino 
carries 5% of the initial proton energy. The flrst break energy in the neutrino spectrum, E^j, 
is due to the break in the photon spectrum. 

£;HL ^ 6.5 X IQS 300 Q^y ^^Q^ 
^br,MeV 

for HL GRBs or E^? = 7.23 x 10^^^^ GeV for LL GRBs. The total energy to be emitted 

^br,keV 

by neutrinos of energy Ei, can be expressed as, 



^•^^^ ^^^^^-^ ^^^^ 



where E^^" is the total isotropic energy of the emitted gamma-ray photons in the energy 
range of IkeV to lOMeV, which is available from the observations [29], ~ 0.3 is the energy 
fraction carried by electrons, and k is a normalization factor considering the differential 
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neutrino spectrum. The relativistic electrons produce the photons by synchrotron radiation 
and inverse Compton scattering of low energy photons, so four orders of magnitude in photon 
energy corresponds to two orders of magnitude in the energy of the radiating charged leptons. 
Photon spectral index 72 = 2.25 corresponds to a spectral index of —2.5 of the relativistic 
electron spectrum. This correspond to k = 1.8 assuming photon fiuence is proportional to 
neutrino luminosity. If the electron spectral index is —2 one would have k — In(lOO) ~ 
4.6. The internal shocks of GRBs are typically collisionless and Coulomb interaction is not 
important between electrons and protons. It is still an open question regarding how some 
energies of Fermi accelerated protons are transferred to electrons and whether electrons and 
protons would have a same spectral index. In principle, the spectral indices of electrons 
and protons in the shocked region could be different. In our calculations, the electron 
spectral index is taken as —2.5 (based on the observed photon spectrum), while the proton 
spectral index is still assumed to be ~ (—2) as has been adopted by all the previous neutrino 
calculations. 

The break (Eq.(lO)) due to the break energy in the photon spectrum is caused by the break 
in (Eq.(6)). If one takes into account muon and pion cooling, a second break appears in 
the neutrino spectrum 

,,dN^{E^) _ 3U (1 - 6e) ^,3, [1 E^<E: 



^^^^ ~ Tli^^^-^ 1 (f)-^ > Ft ^12) 

Pion cooling energy is ten times higher than the muon cooling energy. We would be overes- 
timating the total neutrino flux if we use the pion cooling energy to derive the second break 
energy in the neutrino spectrum. The muon cooling break can be expressed as a function 
of GRB parameters. 

E^'^L ^ 2.56 X lO^e^^e^^/^L-^sf r^oo^.,-3GeV (13) 

for HL GRBs or E^^^ = 3.16 x Wel''^eB^^'^L~%^T\^ty^2GeV for LL GRBs. Here the param- 
eter eB denotes the fraction of internal energy carried by the magnetic fields. We can see 
that if ee and are comparable their contributions are cancelled out. They are assumed 
to be equal to 0.3 in our calculations. Neutrino oscillation only redistributes the number of 
different flavors of neutrinos. The total number of neutrinos remains unchanged. 
The maximum proton energy achievable in HL and LL GRBs can be derived by equating 
the acceleration time scale with the minimum of the dynamic time scale and synchrotron 
cooling time scale of protons [27]. In the comoving frame, the synchrotron cooling time can 
be expressed in terms of proton energy in the source rest frame as, 

atmjEpB^ 

where nip and nie are proton and electron masses respectively, at is the Thompson scattering 
cross section and Be is the comoving magnetic field. The comoving magnetic field strength 

can be expressed as a function of radial distance, luminosity, Lorentz factor the unknown 
shock strength (which depends on the relative Lorentz factor of the two colliding shells), 
and the unknown shock equaiartition parameters [28]. For = 0.3 and reasonable shock 
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strength, one could estimate ^ ~ (2.5 x 10^G)Ll{^ri3 with r^|^ = 0.27r^oot^ _3 for HL and 
B^^ ~ (4.42 X 10^G)L4^Vr3\ rfg^ = 30TIqU^2 for LL GRBs, ris = r/lO^^ jg ^j^g dissipation 
radius where internal shocks are formed. The acceleration time scale in the comoving frame 
in terms of proton energy in the source rest frame is, 



The dynamical comoving time scale is 



T 

tdyn = ■ (16) 



In our calculations, we use ta — m.m{tsync,tdyn) to calculate the maximum proton energy 
in the shock. For typical parameters, one usually has tgync < t^yn, one can then derive the 
expression of the maximum proton energy using the condition tsync = ta, which gives 

^p,max — r OR ' ^^"^^ 

which gives E^^^^ ^ 6 x W^^SM^ GeV for HL GRBs and E^^^^ ^2x 10^°^^ GcV for 

LL GRBs. After averaging over the parameters, including log normal distributions in Lorentz 
factor and variability time, the broken power-law distribution of the luminosity distribution, 
and the redshift distribution of GRBs (see below for detail), we obtain the average values of 
maximum proton energies for LL and HL GRBs. 

The neutrino spectrum emitted from the source is to be corrected for the redshift of the 
GRB to derive the total spectrum to be observed on earth. The total isotropic energy can 
be expressed as a product of the isotropic luminosity (L^) and the duration of the burst (T^). 

The observed neutrino spectrum on earth — ^^^i,*^ ^ is a function of F, L^, T^, z, and E^" 

dNfjEf) _ dN^{E^) 1 

dEi^ - dE^ A^d?{z)^^^'^ ^^^^ 

where r\z) is the comoving radial coordinate distance of the GRB. For a spatially flat universe 
with + Q.jn — 1, the comoving distance can be expressed as 

d{z) = / ^ , ' — (19) 



In our calculations we use Qa = 0.73, flm = 0.27 and the Hubble constant Hq — 
71 km sec~^ Mpc~^ from [30]. When calculating the diffuse neutrino background from GRBs 
we assume that GRBs follow the star formation rate^. The diffuse neutrino flux on earth 
from GRBs distributed up to a redshift of Zmax can be expressed as, 

M?(L„ r, r„ t., E?) ^ J- mm^,n,^ (20) 



"^If however GRBs do not trace the global star-forming rate, our normalization factor needs to be modified. 
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where, the star formation rate is given by [31] 



g3.4. 



Rgrb{z) = 23po ^3.4. + 22 ■ ^^^^ 

The local rates po, of the LL and HL populations have been assumed to be about 
SSOGpc""^^"^ and l.lGpc~^yr~^, respectively [12]. We obtain the neutrino fluxes from 
LL and HL populations distributed up to a distance of Zmax = 5. The neutrino fluxes have 
been averaged with luminosity, Lorentz factor, variability time and burst duration distri- 
bution functions to take into account their variations in the population (the variations in 
luminosity and duration also reflect variation in isotropic energy) 

MK')= T"' M'^\L^,T,Ta,U,E:')S{L^)G{T)P{Ta)Q{U)dL^^^^ 

JLj^i JTi JT^i Jty\ 

Here S{L^) is the broken power law luminosity function (eqn(2)), with Lf^ — 10^^ erg s~^, 
~ 1.15, ~ 2.5 for HL and L^^ = 10^^ erg s-^ a];^ ~ 0.6 and a^^ ~ 4.5 for LL 
GRBs (following [12]). For the HL population the lower and upper limits in luminosities 
are L}^\ = 10^° erg s^^ and = 10^^ erg s~^, respectively, and the corresponding values 
for the LL population are 10"^^ erg s~^ and 5 x 10^^ erg s~^, respectively. The Lorentz factor 
distribution function G(r) is asssumed to be log normal. For HL GRBs, we assume the 
mean of the distribution at log(r^) = 2.6, which correspond to F^ 398, with a standard 
deviation of 0.3. For LL GRBs, we assume the log normal distribution of Lorentz factors is 
peaked at log(Fm) = 1.1 [21] with a standard deviation 0.1. We allow the values of Lorentz 
factors to vary from 100 to 1000 for the HL GRBs, and from 5 to 50 for the LL population. 
The durations of the bursts are also assumed to follow log normal distributions P{Td), with 
the mean at log{Td,m) — 1-47, 3.5; which corresponds to about 30sec and 3000sec, for the 
HL and LL populations, respectively, with standard deviations 0.5 and 0.1. The scattering 
in Amati's relation [23] does not affect our results qualitatively. Since the photon index 
72 does not differ signiflcantly among bursts [18], we take 72 = 2.25 universally for both 
the HL and LL populations. We also consider variation in variability time of GRBs. Their 
distributions are assumed to be log normal, represented by Q{ty) in eqn.(22) with the mean 
at \og{ty^rn) — —1-52, 2.6 and standard deviations 0.3, 0.1 for HL, LL GRBs respectively. The 
maximum energies of neutrinos have also been averaged with the Lorentz factor, luminosity 
and variability time distribution functions as discussed above. 

K%^x = r r E-,rnaxiL-y, F, QS{L^)Gir)Q{QdL^drdt,. (23) 

Also, the redshift correction on the observed neutrino energy has been averaged over the 
entire volume between z = to z = 5. 



j^av,ob ^ r^av RGRB{z)dV {z) / {1 + z) 

v,max J^R^^^(^z)dV{z) ^ ^ 



3 Results and Discussions 

The neutrino spectra from HL and LL populations are plotted in units of GeV cm^^ s^^ sr^^ 
in Fig.l. for 71 = 1, 1.1, 1.3 and 1.5 and 72 = 2.25. We have assumed eg = 0.3 and 63 = 0.3 
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in our calculations. The long dashed (green) lines represent neutrino fluxes from the LL 
population and the solid (red) lines are for the HL population. The Waxman-Bahcall [32] 
limit with z evolution (including all three flavors) has been shown by short dashed (blue) 
lines, and the dotted (pink) lines represent the sensitivity of the IceCube detector [33] after 
three years of operation. The current limit from AMANDA-II data taken between 2000 and 
2003 [34] is also shown by triple dotted hue above the WB upper bound. The HL component 
has an average maximum energy of 6.3 x 10^ GeV and for LL population it is 1.5 x 10^ GcV. 
Above about 10*^ GcV the contribution of LL population to the diffuse neutrino background 
is more than the HL population. The LL component continues to below 100 GeV but 
the HL component has a lower cut-off energy of about 1000 GeV. Between 1000 GeV and 
4 X 10^ GeV their ffuxes are equal. We have also plotted the GZK neutrino ffux from [35] 
in thin dashed lines for a comparison. Comparing with the Icecube sensitivity, we can see 
that the GZK neutrinos are generally not detectable by Icecube, while the contribution of 
the LL-component would be detectable if 71 > 1.1. In this case we expect 30/cm~^sr~^ 
neutrino events per year or more in IceCube at an energy of lO^GeV. If we compare our 
derived neutrino spectra from HL GRBs with some of the previous results [4] the width 
of the intermediate region between the two break energies is narrower. This is because we 
have used the muon cooling energy instead of the pion cooling energy in determining the 
second break energy of the spectra to avoid overestimation of neutrino ffuxes. However, in 
this way we slightly underestimate the ffuxes as they are produced by tt"*" decay unlike Vp, 
and v^. We also ffnd that our results are very sensitive to variations in the mean value of 
Lorentz factor. As the break energies in the neutrino spectra depend on and a slight 
deviation in the value of F changes the neutrino ffuxes significantly. In the future it might 
be possible to set limits on the values of GRB Lorentz factors if the detection or upper limit 
of high energy neutrino background from GRBs is set up. Overall, our results show that the 
previously neglected LL population of GRBs significantly contribute to the diffuse neutrino 
background as their event rate is very high. This conclusion would be strengthened as Swift 
detects more LL GRBs in the years to come. 

Finally we notice that after we submitted the ffrst version of this paper to astro-ph, similar 
results were independently reported by Murase et al.[36]. 

We thank the anonymous referee for constructive reports that help to improve the paper. 
This work was supported by NASA under grants NNG05GB67G and NNG06GH62G. 
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Figure 1: The diffuse neutrino ffuxes expected on earth from low and high luminosity popu- 
lations of GRBs distributed up to a redshift of 5, being plotted in units of GeV cm~^ s^^sr"^ 
against the observed neutrino energy. The four panels are for photon spectral indices 
7i = 1, 1.1, 1.3, 1.5 (clockwise from left bottom), and 72 = 2.25 are adopted for all of them. 
The Waxman-Bahcall limit with z evolution [32] has been shown with short dashed (blue) 
line. The double dotted (pink) line represents IceCube experiment's sensitivity [33] after 
3 years of operation and the current limit from AMANDA-II data (2000-2003) [34] triple 
dotted (orange) line is just above the Waxman-Bahcall limit. These limits include all the 
three flavors of neutrinos. GZK neutrino flux has been plotted in thin dashed line from [35] 
for comparison. 
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